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ABSTRACT 



We present the results obtained by INTEGRAL on 
the Soft-Gamma Ray Repeater SGR 1806-20 . In 
particular we report on the temporal and spectral 
properties of the bursts detected during a moderately 
active period of the source in September and October 
2003 and on the search for quiescent emission. 

Key words: Gamma Rays : bursts; pulsars: general; 
stars: individual (SGR 1806-20). 



1. INTRODUCTION 



Soft Gamma-ray Repeaters (SGRs) are a class of pe- 
culiar high-energy sources discovered through their 
recurrent emission of soft 7-ray bursts. These bursts 
have typical durations of ~0.1 s and luminosit ies in 
the range 10 39 -10 42 ergs s _1 (see lHurlevI {2000) for a 
review of this class of objects). Occasionally, SGRs 
also emit giant bursts that last up to a few hun- 
dred seconds and exhibit rema rkable pulsations tha t 
reveal their spin periods (e.g. iMazets et alJ l)1979[) . 



iHurlev et a.U fiflflflhL 

The bursting activity and the persistent emission 
observed in the ^0.5-10 keV energy range are gen- 
erally explained in t he framework of the "Magne - 
tar" model (see e .g. iDuncan fc Tho mpson! ll!992ft . 
iPaczvhskil l)1992ft . iThompson fc Duncan! l)1995jv Tas 
caused by a highly magnetized (B ~10 15 G) slowly 
rotating (P ~ 5-8 s) neutron star. In this model 
magnetic dissipation causes the neutron star crust 
to fracture. These fractures generate sudden shifts 
in the magnetospheric footpoints, which trigger the 
generation of Alfven pulses, which in turn acceler- 
ate electrons above the pair-production threshold, 
resulting quickly in an optically thick pair-photon 
plasma. The cooling of this plasma generates the 
typical short bursts of soft 7-ray radiation. The 
longer bursts are powered by magnetic reconnection, 
and involve the entire neutron star magnetosphere. 

SGR 1806-20 is one of the most active Soft Gamma- 
ray Repeaters. Here we report new observations 
of this source obtained with the INTEGRAL satel- 
lite in September and October 2003 during a period 
of moderate burstin g activi t y llMereghetti fc Ootzl 
I 20031J: iGotz et all T2q03at IHurlev et all l2003bl: 
iMereghetti et all l2003cHOotz et all l2003bfl . These 
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data have two advantages compared to previous ob- 
servations in the soft 7-ray energy range of bursts 
from this source. First, they have been obtained with 
an imaging instrument, thus we can exclude that the 
bursts originate from a different source in the field. 
Second, they have a good sensitivity and time resolu- 
tion which allows us to study the spectral evolution 
of relatively faint bursts. 



2. BURSTING ACTIVITY 



2.1. Detection and Localization 

In July 2003 a new active peri od of SGR 1806-20 w as 
detected by the IPN (see e.g. iHurlev et all l|2003al) L 
The source remained active throughout the month 
of August and so an INTEGRAL Target of Oppro- 
tunity Observation (TOO) was triggered. INTE- 
GRAL observed it for 240 ks starting on Septem- 
ber 3 2003, and 3 faint bursts were detected. After- 
wards the monitoring of the source continued during 
the Galactic Center Deep Exposure (GDCE) as part 
of INTEGRAL Core Program observations (yield- 
ing an additional exposure of about ~1 Ms on the 
source), during which 21 bursts have been detected. 
All the bursts have been detected and localized in 
near real ti me by the INTEGRAL Bu rst Alert Sys- 
tem (IB AS: IMereghetti et alJ l|2003al) L using IBIS 
l|Ubertini et alll2003|) data in the 15-200 keV range. 
The burst detection rate vs. time is shown in Fig. ^ 
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Figure 1. Burst detection rate per day. The shaded 
areas represent the time periods when the source 
was outside the IBIS field of view (short gaps due 
to perigee passages and satellite slews are not indi- 
cated). 

All the bursts have been localized within 2 arcmin- 
utes fro m the X-ray positio n of the quiescent coun- 
terpart <|Kanlan et all 12002) . Sinc e the typical error 
circle is 2.5 arcmin (90% c.L, see IMereghetti et al.l 
(2004)), we are confident that all the bursts originate 
from SGR 1806-20 and not from the possible newly 
discovered Soft Gamma-Ray Repeater SGR 1808-20 



l|Lamb et all l2003[) . which is located at 15 arcmin 
from the source. 



2.2. Data Analysis 

We have analyzed only ISGRI (iLebrun et all l2003|) 
data, since PICsIT l|Labanti et all 2003) in its stan- 
dard operation mode does not have the combined 
time resolution and sensitivity to study such short 
bursts. The background subtracted IBIS/ISGRI 
light curves of the bursts, binned at 10 ms, are shown 
in Fig. 13 In order to increase the signal-to-noise- 
ratio, they were extracted from ISGRI pixels illumi- 
nated by the source for at least half of their surface 
and selecting counts in the 15-100 keV energy range 
(most of the bursts had little or no signal at higher 
energy). 

The bursts were detected at various off-axis angles, 
ranging from 2.5 to 13.3 degrees, corresponding to 
a variation of 80% in the instrument effective area. 
The light curves shown in Fig. [21 have been cor- 
rected for this vignetting effect. The total num- 
ber of net counts actually recorded for each burst 
is indicated in the corresponding panel. The burst 
positions in instrumental coordinates are plotted in 
Fig- 01 (see also Tab. H). As can be seen most of 
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Figure 3. Burst positions in instrumental coordi- 
nates. The fields of view of IBIS (large square), SPI 
(hexagon), JEM-X (circle) and OMC (small square) 
are plotted. The units of the grid are degrees. 

the bursts are located in the partially coded field of 
view of IBIS. A fe w bursts fall within the JEM-X 
l|Lund et all |2Q03) field of view, but unfortunately 
they are the faintest ones. The only burst for which 
a simultaneous IBIS/ JEM-X detection has been ob- 
tained (see Fig. 0}, is the one with instrumental 
coordinates Y=-0.97° and Z=-2.22° (burst number 
8). This burst has a rather small fluence: 2.2 xlO -8 
erg cm" 2 (15-100 keV). 
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Burst num. 1 



-54150:53:43 - Tot. Cts: 109 



Burst num. 2 



04T01:15:56 - Tot. Cts: 122 




Burst num. 3 ■ 



-05T07:46:1 5 - Tot. Cts: 131 




Burst num. 4 - 1 0-06T1 8:51 :27 - Tot. Cts 




Burst num. 5 - 10-08T03:27:C3 - Tot. Cts: 770 



Burst num. 6 - '0-08T03:27:56 - Tot. Cts: 562 





Burst num. 7 - 10-08T20:22:06 - Tot. Cts: 250 



Burst num. 8 - 10-1 1T01 :26:27 - Tot. Cts: 147 



Burst num. 9 - 10-14T14:31:45 - Tot. Cts: 302 






Burst num. 10 - 10-15T02:35:19 - Tot. Cts: 162 




Burst num. 11 - 10-15T04:55:00 - Tot. Cts: 777 




Burst num. 12 - 10-15T07:37:47 - Tot. Cts: 210 




Burst num. 13 - 10-15T07:55:53 - Tot. Cts: 203 




Burst num. 14 - 10-15T07:56:21 - Tot. Cts: 259 




Burst num. 15 - 10-15T10:19:31 - To;. C:s: 1"4 



Burst num. 16 - 10-15T1 1:00:28 - Tot. Cts: 131 




Burst num. 17 - 10-15T1 1:17:05 - Tot. Cts: 124 



Burst num. 18 - 10-15T1 2:08:26 - Tot. Cts: 195 



Burst num. 19 - 10-1 5T1 2:37:32 - Tot. Cts: 151 



Burst num. 20 - 10-15T12:51:33 - Tot. Cts: 193 



Burst num. 21 - 1 0- 1 5T1 3: 1 8:26 - Tot. Cts: 409 




40 - 

20^ 

-28 




Burst num. 22 - 1 0-15715:55:1 2 - Tot. Cts: 1145 



Burst num. 23 - 10-15T19:23:38 - Tot. Cts: 561 



Burst num. 24 - 10-15T21 :47:07 - Tot. Cts: 212 






Figure 2. IBIS/ISGRI background subtracted light curves of the SGR 1806-20 bursts in the 15-100 keV range. 
Each panel corresponds to a time interval of one second and the time bins are of 10 ms. Units of the axes are 
time in seconds and vignetting- corrected counts per bin. Time—0 corresponds to the the starting time of the Tqq 
computation and is reported on top of each panel (month - day - UT time) together with the total number of net 
counts. 



None of the bu rsts has been found in a preliminary 
analysis of SPI ijVedrenne et al.ll2003|) data, but fur- 
ther analy sis is underway. Only one burst fell into 
the OMC l|Mas-Hesse et all l2003|) field of view, but 
the large extinction towards the object (A^ ~ 30 
mag) pr events the detection of any optical emission 
(see also lCastro-Tirado et alJ l)1998t) V 



The light curves shown in Fig. have shapes typ- 
ical for SGR bursts. From the light curves we de- 
termined the Tqq duration of each burst (i.e. the 
time during which 90% of the total burst counts are 
accumulated). The Tgo values, reported in Tab. ^ 
range typically from ^0.1 to ^0.2 s for single peaked 
bursts and can be as long as ~0.6 s for double peaked 



4 



Burst 


Fluence 


Peak Flux 




Net 


Off- Axis Angle Y 


Off Axis- Angle Z 


Number 


10~ 8 erg cm~ 2 


10~ 7 erg cm -2 s _1 


ms 


Counts 


degrees 


degrees 


1 


1.65 ± 0.16 


3.74 ± 0.75 


437 


109 


-3.89 


-1.86 


2 


1.84 ± 0.17 


4.26 ± 0.80 


69 


122 


-3.96 


0.15 


3 


1.97 ± 0.17 


2.15 ± 0.57 


579 


131 


4.05 


0.17 


4 


4.92 ± 0.27 


14.77 ± 1.49 


109 


85 


10.26 


8.21 


5 


65.35 ± 0.99 


60.73 ± 3.02 


179 


770 


3.36 


-12.96 


6 


46.58 ± 0.84 


50.42 ± 2.75 


199 


562 


3.35 


-12.92 


7 


5.49 ± 0.28 


15.73 ± 1.53 


99 


250 


3.64 


-7.56 


8 


2.21 ± 0.18 


3.67 ± 0.74 


169 


147 


-0.97 


-2.22 


9 


6.15 ± 0.30 


12.71 ± 1.38 


89 


302 


-6.32 


-4.49 


10 


2.89 ± 0.21 


5.27 ± 0.89 


139 


162 


-5.61 


0.77 


11 


20.99 ± 0.56 


36.40 ± 2.34 


169 


777 


-0.71 


8.87 


12 


3.17 ± 0.22 


9.18 ± 1.17 


129 


210 


-3.69 


3.77 


13 


3.06 ± 0.21 


4.83 ± 0.85 


89 


203 


-3.69 


3.73 


14 


3.89 ± 0.24 


4.86 ± 0.85 


489 


259 


-3.69 


3.78 


15 


3.94 ± 0.24 


15.62 ± 1.53 


69 


114 


-10.03 


-3.26 


16 


3.50 ± 0.23 


7.95 ± 1.09 


99 


131 


-8.63 


-1.40 


17 


2.72 ± 0.20 


7.21 ± 1.04 


189 


124 


-7.29 


0.50 


18 


3.62 ± 0.23 


9.45 ± 1.19 


89 


195 


-5.95 


2.44 


19 


2.50 ± 0.19 


8.49 ± 1.13 


119 


151 


-4.66 


4.42 


20 


3.71 ± 0.24 


14.33 ± 1.47 


159 


193 


-3.45 


6.41 


21 


7.88 ± 0.34 


8.15 ± 1.11 


475 


409 


-3.46 


6.42 


22 


27.21 ± 0.64 


43.71 ± 2.56 


188 


1145 


-3.87 


8.10 


23 


25.17 ± 0.61 


54.80 ± 2.87 


288 


561 


-11.12 


-2.48 


24 


4.29 ± 0.25 


4.24 ± 0.80 


159 


212 


-5.67 


5.12 



Table 1. Fluences (15-100 keV), Peak Fluxes (over 10 ms, 15-100 keV), Tqq durations, net counts (background 
subtracted, but not corrected for vignetting) and instrumental coordinates for all the bursts. 



bursts. In fact the Tqq values of these bursts include 
the "interpulse" period. Some bursts are preceded 
by a small precursor. 

The peak flux and fluence for each burst were first 
derived in counts units from the light curves of 
Fig|21 and then converted to physical units adopt- 
ing a constant conversion factor derived from the 
spectral analysis of the brightest bursts (see below). 
The resulting 15-100 keV peak fluxes and fluences 
are reported in Tab. ^ an d are respectively in the 
range (4-50) xlO -7 erg cm~ 2 s _1 (Ai=10 ms) and 
(2-60) xlO -8 erg cm -2 . Their integral distributions 
are shown in Fig. |SJ Within the large uncertain- 
ties, the fluence distributio n is consistent wi th the 
power law slope found bv lGogiis et alJ l)2000[) using 
RXTE data. Many of these bursts are the faintest 
ever imaged from SGRs at these energies. 

For the bursts with more than 500 net counts wc 
could perform a detailed spectral analysis. The 15- 
200 keV spectra, integrated over the whole duration 
of each burst, were well fitted by an Optically Thin 
Thermal Bremsstrahlung (OTTB) model, which is 
the spectral model that typically provides good fits 
to SGR bursts above 20 keV. Using the latest avail- 
able response matrices, we measured temperatures 
(kT) between 32 to 42 keV. These values are in the 
typical r ange of tempe r ature s found in the literature 
(see e.g. lAtteia et all (^987)) for short bursts from 
SGR 1806-20 . We tried other models, like a power 
law or a black body, but they were clearly ruled out. 
Small known calibration uncertainties at low energies 



(< 40 keV) do not affect much our spectra, whose er- 
rors are dominated rather by statistical uncertainties 
than by systematic ones. The OTTB fit to the spec- 
trum of burst number 22 is shown as an example in 
Fig. El 

Adopting a temperature fcT=38 keV (consistent with 
the average spectra of the brightest bursts) we de- 
rived a conversion factor of 1 count s _1 = 1.5x 10~ 10 
erg cm~ 2 s _1 (15-100 keV), which we adopted for all 
the bursts. 

iGotz et all l)2004ah have used these data, to investi- 
gate the time evolution of the bursts spectra, using 
time-resolved hardness ratios (H R = (H — S)/(H + 
S), where H and S are the background subtracted 
counts in the ranges 40-100 keV and 15-40 keV re- 
spectively). Analyzing the 12 bursts with more than 
200 net counts, they found that some bursts show a 
significant spectral evolution, while others, particu- 
larly those with a "flat topped" profile, do not. Some 
examples are given in Fig. In addition the vari- 
ation of the hardness ratio versus intensity (/) has 
been investigated. Cons idering all the time b ins of all 
the bursts (see Fig. ISjl. lGotz et alJ ()2004a[) found a 
hardness-intensity anti-correlation. The linear corre- 
lation coefficient of the data corresponds to a chance 
probalility smaller than 10~ 3 of being due to uncor- 
rected data. In addition, according to an F-test, 
the data are significantly (at a 5.2 a level) better 
described by a linear fit than by a constant value. 

We have verified that there is no correlation between 
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Burst § 11 Burst # 22 Burst # 5 Burst § 9 




Time [s] Time [s] Time [s] Time [s] 



Figure 7. 15-40 keV light curve (Top PanelsJ, 40-100 keV light curve ('Middle Panels,), and time resolved 
hardness ratio ('Bottom Panels,) for four bursts with good statistics. The time resolved ha rdness ratio for bursts 
number 11,22,9 is inconsistent with a constant value at ~3.5 a level \G6tz et al\ . \200~I\) . 



the hardness ratios of the individual bursts and the 
off-axis angle at which they have been detected. In 
adition we have verifie d that the vignetti ng: correc- 
tion procedure used in lOotz et alJ l|2004afl is consis- 
tent with the values of the flux of the Crab Nebula 
measured in different positions in the field of view of 
IBIS. 



3. PERSISTENT EMISSION 



Persistent (quiescent) emission from SGR 1806- 
20 has been discq r yered at X-ray (<10 keV) energies 
l|Murakami et all ^994). Up to now no detection 
at higher energies has been reported for any of the 
SGRs. We have analyzed IBIS, JEM-X and SPI data 
of the TOO (240 ksec) , but we did not find convinc- 
ing evidence of quiescent emission. 

On the other hand analyzing IBIS/ISGRI Core Pro- 
gram data (~1 Msec of exposure on SGR 1806-20 ), 
we detect the source at ^6 a level in IBIS /ISGRI (see 
Fig. EJ) in the 20-40 keV band. The flux is -3 mCrab, 
consistent with an extrapolation of the spectrum 



measured at lower energies l)Mereghetti et all [2000) . 
Although a detection with such a relatively small 
significance is also compatible with background sys- 
tematic noise in the significance maps, the positional 
coincidence with the X-ray counterpart strengthens 
the detection. 

Also SPI data have been searched for quiescent 
emission. 482 pointings around the position of 
SGR 1806-20 (Core Program and TOO data) have 
been used, yieldin g a total exposure of ~0.798 Ms. 
For the SPIROS (jSkmner k ConnellL l200l analy- 
sis tool the catalogue of the sources detected by 
IBIS/ISGRI has been used as an input. A sig- 
nificance map was generated in the energy range 
between 28 and 48 keV. The measured flux is 
0.69±0.15xl0" 3 photons cm~ 2 s" 1 . This flux (~ 
10 mCrab) is slighly higher than the one measured 
by IBIS, but considering that the SPI detection level 
is just ~ 4.4 a, and that systematic errors are still 
present at this stage in the analysis of both instru- 
ments, we can consider the two results as consistent 
with each other. 

Further GCDE observations during INTEGRAT 
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JEM X 




ISGRI 40-'C0 keV PIF>0.5 




Figure 4- Light curves of the SGR burst detected also 
by JEM-X. 



AO-2 will surely help to further assess the detection 
and the source flux. 

Concerning JEM-X Core Program data, the analysis 
has still to be completed and will surely be helpful 
in order to better asses the possible IBIS and SPI 
detection. 

The detection of pulsed emission is at the moment 
below our sensitivity. 



4. CONCLUSIONS 



Figure 5. Integral distributions of Peak flux (At=10 
ms, Upper Panel,) and fluences ('Lower Panel/ Both 
refer to the 15-100 keV energy range. 



SGR 1606-20: burst § 22 




We can summarize our results as follows: 



• For the first time good evidence for spectral 
evolution of weak SGR bursts is presented and 
a hardness-intensisty a nti-correlation with in the 
bursts has been found ijGotz et alll20044 - This 
results represent a new challenge for the Magne- 
tar model, which predicts that the effetive tem- 
perature of the burst should vary weakly during 
the bursts, while we detect large spectral vari- 
ations whitin the bursts as for number 9 and 
22. 

• We report for the first time the possible detec- 
tion, in IBIS and SPI data, of the quiescent 
emission above 10 keV. 

• IBIS is a very sensitive detector for SGR bursts, 
with fluences down to ^10~ 8 erg cm -2 in imag- 



20 50 100 

channel energy (keV) 

Figure 6. Time averaged spectrum of burst # 22 
(IBIS/ISGRI). 



ing mode. In fact in just 1.5 months it detected 
twice as many bursts from this source as BATSE 
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l2004bf) and the ~0.5 s long faint burst detected at 
13:02:06.33 UT on March 16 (IBAS Alert n. 1677). 
The alert messages have been distributed just —20 
seconds after the bursts, allowing for prompt follow- 
up observations at other wavelenghts. This feature 
will be particularly useful in case of major bursts like 
the ones from SGR 1900+14 or SGR 0526-66. 
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